function [adjMatrix,nodeNames] = getDAG()
nodeNames{1} = 'N5';
nodeNames{2} = 'N9';
nodeNames{3} = 'N19';
nodeNames{4} = 'N20';
nodeNames{5} = 'N27';
nodeNames{6} = 'N29';
nodeNames{7} = 'N34';
nodeNames{8} = 'N38';
nodeNames{9} = 'N26';
nodeNames{10} = 'N28';
nodeNames{11} = 'N40';
nodeNames{12} = 'N41';
nodeNames{13} = 'N57';
nodeNames{14} = 'N42';
nodeNames{15} = 'N63';
nodeNames{16} = 'N64';
nodeNames{17} = 'N65';
nodeNames{18} = 'N66';
nodeNames{19} = 'N67';
nodeNames{20} = 'N68';
nodeNames{21} = 'N33';
nodeNames{22} = 'N35';
nodeNames{23} = 'N37';
nodeNames{24} = 'N39';
nodeNames{25} = 'N43';
nodeNames{26} = 'N44';
nodeNames{27} = 'N69';
nodeNames{28} = 'N70';
nodeNames{29} = 'N71';
nodeNames{30} = 'N73';
nodeNames{31} = 'N51';
nodeNames{32} = 'N52';
nodeNames{33} = 'N88';
nodeNames{34} = 'N89';
nodeNames{35} = 'N90';
nodeNames{36} = 'N91';
nodeNames{37} = 'N59';
nodeNames{38} = 'N141';
nodeNames{39} = 'N143';
nodeNames{40} = 'N146';
nodeNames{41} = 'N147';
nodeNames{42} = 'N150';
nodeNames{43} = 'N152';
nodeNames{44} = 'N153';
nodeNames{45} = 'N187';
nodeNames{46} = 'N81';
nodeNames{47} = 'N84';
nodeNames{48} = 'N113';
nodeNames{49} = 'N114';
nodeNames{50} = 'N115';
nodeNames{51} = 'N116';
nodeNames{52} = 'N118';
nodeNames{53} = 'N119';
nodeNames{54} = 'N120';
nodeNames{55} = 'N121';
nodeNames{56} = 'N122';
nodeNames{57} = 'N100';
nodeNames{58} = 'N104';
nodeNames{59} = 'N108';
nodeNames{60} = 'N112';
nodeNames{61} = 'N117';
adjMatrix = zeros(length(nodeNames));
adjMatrix(39,1) = 1;
adjMatrix(45,1) = 1;
adjMatrix(1,2) = 1;
adjMatrix(3,2) = 1;
adjMatrix(27,2) = 1;
adjMatrix(5,4) = 1;
adjMatrix(6,4) = 1;
adjMatrix(27,4) = 1;
adjMatrix(31,4) = 1;
adjMatrix(7,8) = 1;
adjMatrix(3,9) = 1;
adjMatrix(6,10) = 1;
adjMatrix(1,11) = 1;
adjMatrix(3,11) = 1;
adjMatrix(5,12) = 1;
adjMatrix(6,12) = 1;
adjMatrix(7,13) = 1;
adjMatrix(1,14) = 1;
adjMatrix(3,14) = 1;
adjMatrix(5,14) = 1;
adjMatrix(6,14) = 1;
adjMatrix(31,14) = 1;
adjMatrix(1,15) = 1;
adjMatrix(3,15) = 1;
adjMatrix(27,15) = 1;
adjMatrix(31,15) = 1;
adjMatrix(5,16) = 1;
adjMatrix(6,16) = 1;
adjMatrix(1,17) = 1;
adjMatrix(3,17) = 1;
adjMatrix(3,18) = 1;
adjMatrix(6,18) = 1;
adjMatrix(27,18) = 1;
adjMatrix(5,19) = 1;
adjMatrix(6,19) = 1;
adjMatrix(27,19) = 1;
adjMatrix(31,19) = 1;
adjMatrix(7,20) = 1;
adjMatrix(1,21) = 1;
adjMatrix(1,22) = 1;
adjMatrix(5,23) = 1;
adjMatrix(3,24) = 1;
adjMatrix(6,25) = 1;
adjMatrix(7,26) = 1;
adjMatrix(28,31) = 1;
adjMatrix(30,31) = 1;
adjMatrix(31,32) = 1;
adjMatrix(28,33) = 1;
adjMatrix(27,34) = 1;
adjMatrix(27,35) = 1;
adjMatrix(27,36) = 1;
adjMatrix(27,37) = 1;
adjMatrix(29,37) = 1;
adjMatrix(39,38) = 1;
adjMatrix(40,38) = 1;
adjMatrix(41,38) = 1;
adjMatrix(43,38) = 1;
adjMatrix(44,38) = 1;
adjMatrix(42,45) = 1;
adjMatrix(38,46) = 1;
adjMatrix(38,47) = 1;
adjMatrix(38,48) = 1;
adjMatrix(27,49) = 1;
adjMatrix(31,50) = 1;
adjMatrix(29,51) = 1;
adjMatrix(1,52) = 1;
adjMatrix(3,53) = 1;
adjMatrix(5,54) = 1;
adjMatrix(6,55) = 1;
adjMatrix(7,56) = 1;
adjMatrix(27,57) = 1;
adjMatrix(28,58) = 1;
adjMatrix(27,59) = 1;
adjMatrix(30,60) = 1;
P = topologicalPermutation(adjMatrix);
adjMatrix = adjMatrix(P,P);
nodeNames = nodeNames(P);
